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Abstract 

On the basis of the approximate dynamical equations describing the behavior of quark- 
gluon plasma (QGP) in the semiclassical limit and Yang-Mills equation, the kinetic equa- 
tion for longitudinal waves (plasmons) is obtained. With the Ward identities the gauge 
invariance of obtained nonlinear Landau damping rate is proved. The physical mecha- 
nisms defining nonlinear scattering of a plasmon by QGP particles are analyzed. The 
problem on a connection of nonlinear Landau damping rate of longitudinal oscillations 
with damping rate, obtained in the framework of hard thermal loops approximation, is 
considered. It is shown that the gauge-dependent part of nonlinear Landau damping rate 
for the plasmons with zero momentum vanishes on mass-shell. 
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1. INTRODUCTION 



In recent 15-20 years, a theoretical investigations of properties of quark-gluon plasma 
has been of great interest. It is connected with intensive looking for a QGP in the 
experiments with collision of ultrarelativistic heavy ions and application concerning the 
physics of the early universe. 

Two methods to study of the nonequilibrium phenomena in a QGP are used: method of 

temperature Green functions and kinetic approach. Significant progress has been achieved 
in the development of the first method. The effective perturbative theory was constructed 
in the papers [1] by Pisarski and Braaten, Frenkel and Taylor on the basis of the resum- 
mation of so-called hard thermal loops (HTL's), and the problem on the sign and gauge 
dependence of the damping rate of the long wavelength excitations in QGP was solved. In- 
dependent solution of "the plasmon puzzle" problem was received by Kobes, Kunstatter, 
and Rebhan [2] . The damping rate for heavy quarks interacting with light thermal quarks 
and gluons, was found in [3], for soft plasmon and plasmino - in [4, 5]. The damping rate 
for energetic fermions and bosons was derived by Lebedev and Smilga [6] with extension 
of the program of resummation outlined in Ref . [1] , by means of inclusion of higher order 
effects in a hard propagator. The alternative derivation of damping rate for energetic 
fermions, based on introduction of an infrared cutoff was deduced by Burgess, Marini, 
and Rebhan [7]. The progress of the thermal QCD makes possible a new look at the 
existing of kinetic theory of QGP, developed by Heinz, Winter, Elze, Vasak and GyuUasy, 
Mrowczyhski and others [8], and it has given impetus to its further development. 

In spite of the fact that the language and methods of these aproaches are very different, 
there are close similarities between HTL approach and transport theory. Originally the 
kinetic theory was used by Sihn to derive HTL in the photon's self-energy [9]. HTL's 
in the quark and gluon self-energies can be computed similarly. Moreover, Kelly, Liu, 
Lucchcsi and Manuel [10] have shown that the generating functional of HTL's (with an 
arbitrary number of soft external bosonic legs) can be derived from the classical kinetic 
theory of QGP. This points to the classical nature of the hard thermal effects. For hard 
excitations, the damping rate has been computed by Heiselberg and Pethick [11] from a 
Boltzmann - like equation, with the collision terms included. 

A further step in development of kinetic theory was made by Blaizot and lancu [12] . In 

contrast to the early papers on transport theory of QGP [8], these authors use from outset 
the ideas developed in thermal QCD in deriving of the kinetic equations. The equations 
obtained by them on the basis of a truncation of the Schwinger-Dyson hierarchy isolate 
consistently the dominant terms in the coupling constant 51 in a set of equations, which 
describe the response of plasma to weak and slowly varying disturbances, and encompass 
all HTL's (conserning recent investigations in this area, see also papers [13-16]). However, 
here it should be noted that if the influence of the average fermionic field is neglected, 
then the expression for current induced by soft gauge fields, obtained in [12] (and the 
nonlinear equation of motion, connected with it) fully coincide with the corresponding 
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expression obtained in [10] from usual classical kinetic theory on the basis of consistent 
expansion of distribution function in powers of the coupling constant. This somewhat 
justifies use of the (semi)classical kinetic equations found in [8], in spite of the fact that 
intermediate approximation schemes in which these equations were derived, mix leading 
and nonleading contributions with respect to the powers of g and so, are not entirely 
consistent. 

Such close interlacing of two methods of investigation of nonequilibrium phenomena 
in QGP leads to the question: can we calculate the damping rate of soft bosonic modes 
corresponding to the hard thermal result [4], remaining in the framework of classical 
(semiclassical) kinetic theory only? Xiaofei and Jiarong [17] were the first to put this 
question. Because of obtained results, they have given a positive answer. 

As was shown by Heinz and Siemens [18], in linear approximation the Landau damping 
is absent in QGP. In fact, the only mechanism, with that one can connect the damping 
following from the kinetic theory with one from HTL approach, is so-called nonlinear 
Landau damping. It bounds up with the nonlinear effects of waves interaction and par- 
ticles in QGP. The multiple time-scale method which has proved successful in study of 
the nonhnear properties of electromagnetic (Abehan) plasma [19], was used in [17] for 
determination of this association. By means of this method the nonlinear shift of the 
mass-shell of the longitudinal modes in the temporal gauge has been obtained by Xiaofei 
and Jiarong. Its imaginary part defines required nonlinear Landau damping rate. Futher, 
the limiting expression of the derived damping rate for k = 0-mode was obtained, and 
numerical computations for approximate estimate were performed. The value derived by 
this means is in close agreement with similar numerical one obtained by Braaten and 
Pisarski [4] in the framework of effective perturbative theory. 

However, under close examination of above-mentioned paper we found certain mistakes 
in computations, which were of both principle and nonprinciple character. As it was shown 
in our early paper [20], the elimination of these inaccuracies finally leads not only to a 
numerical modification of the limiting value of nonlinear Landau damping rate obtained in 
[17], but what is more important, it changes the sign of obtained expression. This points 
to some prematurity of the statements in [17] on obtained connection between nonlinear 
Landau damping rate and damping rate from the HTL-approach. 

In this paper, that is further development the ideas outlined in Ref . [20] , we consider 
the above problem, using the approach based on obtaining of kinetic equation for waves 
in quark-gluon plasma developed by Kadomtsev, Silin, Tsytovich and others [21] in con- 
nection with ordinary plasma. We have shown that nonlinear Landau damping rate 7'(k) 
(which even is not of fixed sign for arbitrary value of k) for longitudinal waves in QGP 
defines two various processes: the effective spectral pumping of energy from short to long 
waves (with complete conservation of excitations energy), and properly nonlinear dissi- 
pation of plasma waves energy in the medium. The main conclusion of this work is that 
there is a need to compare the piece of 7'(k) that corresponds to nonlinear dissipation 
of waves energy and is a positive for any value of a wave vector k and, in particular for 
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k = - mode, with damping rate of soft boson modes from HTL-approximation. 

The outhne of the paper is as follows. In Sec. 2 we derive a system of self-consistent 
equations in the covariant gauge for regular (coherent) and random parts of both the 
distribution functions of QGP particles and gauge field. In Sec. 3 the first order ap- 
proximation of the colour current is considered and the correlation function of random 
oscillations is introduced. In Sec. 4 we discuss the consistency with gauge symmetry of 
used approximation scheme. In sec. 5 the second and the third orders approximation of 
the colour current are studied, and the terms leading in the coupling constant are sepa- 
rated. In Sec. 6 the kinetic equation for longitudinal waves in QGP is derived. In Sec. 
7 the nonlinear Landau damping rate is rewritten in the term of HTL-amplitudes. In 
Sec. 8 by the effective Ward identities, the gauge invariance of obtained damping rate is 
proved. In Sec. 9 the physical mechanisms defining the nonlinear scattering of waves by 
plasma particles are considered. In Sec. 10 association of the nonlinear Landau damping 
rate with damping rate, obtained on the basis of hard thermal loops approximation is 
discussed. In Sec. 11 the estimation of a value of the nonlinear Landau damping rate of 
plasmons at rest (at vanishing three-momentum) is made. In Sec. 12 it is shown that 
the gauge-dependent piece of nonlinear Landau damping rate vanishes on mass-shell. In 
Conclusion possible ways of further development of the scrutinized theory are discussed. 



2. THE INITIAL EQUATIONS. THE RANDOM-PHASE APPROXIMA- 
TION 



We use metric g^'^ = diag{l, —1, —1, —1) and choose units such that c = ks = 1. The 
gauge field potentials are Nc x A'^c-iiiatriccs in a color space defined by = A^t" with 
A^^ — 1 hermitian generators of SU{Nc) group in the fundamental representation. The 
field strength tensor = F^^t"' with 

f;, = d,Al - d^Al + gr'-A'^Al (2.1) 
obeys the Yang-Mills (YM) equation in a covariant gauge 

d,F^\X) - ig[A,{X), F^'^iX)] - C'd'^d^A^iX) = -f{X), (2.2) 
where ^ is a. gauge parameter, is the colour current 

/ = ge J d'pp-'iSpt^f, - /,) + Tr (r%)], (2.3) 

where are hermitian generators of SU(Nc) in the adjoint representation ((T")'"^ = 
—if°'^,Tr{T"T^) = NcS"'''). Wc denote the trace over color indices in adjoint representa- 
tion as Tr. Distribution functions of quarks fg, antiquarks fg, and gluons fg satisfy the 
dynamical equations which in the semiclassical limit (neglecting spin effects) are [8] 

1 df - 

p'-^.Ua ± l^9J^{F,., = 0, (2.4) 
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1 df 
I^V.f, + -9P^{T,., ^} = 0, 

where and are covariant derivatives which act as 

V,^d,-ig[A^{X),-l 

V^ = d^-ig[A^{X),-l 

[ , ] denotes commutator, {, } denotes anticommutator, and T^i, are defined as An — 
A'^^T"', T^v = F^JT"-. Upper sign in the first equation (2.4) refers to quarks and lower one 
- to antiquarks. 

We begin with consideration of dynamical equations (2.4). The distribution functions 
fq^q and fg can be decomposed into two parts: regular and random ones, where latter are 
generated by spontaneous fiuctuations in the plasma 

/. = /f + /J, s = qA.9. (2.5) 

so that 

(/.) = , (/J) = 0. (2.6) 

Here, angular brackets (■) indicate a statistical ensemble of averaging. The initial values of 
parameters which characterize the collective degree of a plasma freedom is such statistical 
ensemble. For almost linear collective motion to be considered below it may be initial 
values of oscillation phases. 

Further we set 

A, = A^^+Al, {AD = 0, (2.7) 

by definition. The regular (background) part of the field A^ will be considered equal to 
zero. The condition for which the last assumption holds, will be closer considered in Sec. 4. 

Averaging the equation (2.4) over statistical ensemble, in view of (2.5)-(2.7), we obtain 
the equations for the regular parts of the distribution functions f^^ and 

^-^9P'{{{FI)nl.^]). (2.8) 



P'd.f^ = ^9P'{[Al, /J]) - 19p'{{{:fI)l, ^}) - Igpna^^NL), ^} 



4.P''({(^J.W,g}). 
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Here, indices "L" and "NL" denote the linear and nonlinear parts with respect to field 
of the strength tensor (2.1). 

Subtracting (2.8) from (2.4), we define the equations for fj^ and /J 



^\9T^{{{FI)nl. g} - ({(FjJiVL, g})), (2.9) 

P'd,f^ - igp'iiAl, /J] - {[Al, /J])) - l9P^{{:Fp„ 
1 dF dF 1 

-2^P^({(^J.)l, ^} - ({(^i)L, ^})) - :^gp'm.)NL - ((^JjiVL), 

1 df^ 

—^9P'{m^)NL, ^} - {m.)NL, ^}». 

The system of equations (2.8) and (2.9) is suitable for investigation of nonequilibrium 
processes in QGP such that the excitation energy of waves is small quantity in relation to 
the total energy of particles. In this case it can be used expansion in powers of oscillations 
amplitude of the random functions /J 

oo 

/J=E/J^"^ ^=^'5^'^' (2-10) 

n=l 

where /J*^"-' collects the contributions of the n-th power in A^. Substituting expansion 
(2.10) into (2.9), and collecting terms of the same order in ^4^, we derive the system of 
equations 

p'^f.^' =T\9P'm.)J-^}^ (2.11) 

P'd.flf^ = i9jf{[Al flf\ - ([AJ, /J«]))T (2.12) 
Ti^9P'{{{FI)l. ^} - ({(Fjj^, ^})) T '-9P'{{Fl)^L - {{F^M, -^h 
P'djjf = ^9P^([Al, f^f] - {[Al, /ff )])) T 19P^{{{FP„ (2.13) 

-m.)L, ^})) T ^9P'{m.)NL, ^} - {{{Fpr^L, ^})) etc.. 
Similar equations are obtained for /J^*-"^ , n = 1, 2, 3, . . . . 
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The nonlinear colour current is expressed as 

oo 

j.-j^+jj:, (j.)=tf ,^=E^^'^^ (2-14) 

n=l 

where 

iJW = gf^ J Ap.[Spt"(/J(") - Z^"^) + IV (TVJ^"^)]- (2.15) 

Now we turn to the Yang-Mills equation (2.2), connecting the gauge field with the 
colour current. Averaging Eq. (2.2) and subtracting the averaged equation from (2.2) in 
view of Eqs. (2.7) and (2.14), we find (for = 0) 

d^,{F^'"')L - r^d^d^Al + f^^^" = (2.16) 

+i9{[Al, {F^nL] - {[Al, [F^nL])) + 9\[Al, A^l] - {[Al, [A^^, A^l])). 

Here, in the l.h.s. we collect all linear terms with respect to Aj^ and we denote: j^'^ = 

jT(2)u jTi3)u _|_ account for nonfinear interaction between waves and particles in 

QGP (in first non- vanishing approximation over the energy of waves) , it is sufficiently to 
restrict the consideration to third order in powers of in expansion (2.10). 

We introduce the following assumption. Eqs. (2.8) represent the kinetic equations for 
averaged distribution functions. The correlation functions in the r.h.s. of these equations 
have meaning of the collision terms due to particlc-wavc interaction and describe the in- 
fiuence of plasma waves to a background state. Recently research of similar equations has 
attracted detailed attention [16], since on scale large wavelengths of collective excitations 
(A ~ l/g^T) they lead to Bodeker's effective theory [13]. We suppose that a characteristic 
time of nonlinear relaxation of the oscillations is small quantity as compared with a time 
of relaxation of the distribution particles f^. Therefore we neglect by change of regular 
part of the distribution functions with space and time, assuming that these functions are 
specified and describe the global equihbrium in QGP 

fR _ 2iV;g(po) ,N I fR= fO^ o 2g(po) ,N I (2 17) 

Jq,q — Jq,q ^ (^27:^ ^ g(p«)/TTM + l ' •'s ~ a (^27r)^ e^^")/"^ - 1 ' ^ ' 

where Nf - being the number of flavours for massless quarks, is the four-velocity of 
the plasma at temperature T, and 11 is the quark chemical potential. 

3. THE LINEAR APPROXIMATION. THE CORRELATION FUNCTION 
OF THE RANDOM OSCILLATIONS 

We will now come to the derivation of kinetic equation for waves. The initial equation 
is Eq. (2.16). The l.h.s. of Eq. (2.16) contains a linear approximation of the colour 
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current, explicit form of which is easily defined from Eq. (2.11). We prefere to work in 
momentum space; the corresponding equations are obtained by using 



ikX 



and similar translations for fj^, /J. The result of Fourier transformation for Eq. (2.11) 
is 

fq,q {k,p) = Tg-j——. — -^A^{k) , >{k,p) = -g-j——. — ■^A^ik), (3.1 
pk + zpoe op^ pk + zpoe op^ 

e +0. 

Here x''^{k,p) = {pk)g''^ -p^k^. Substituting (3.1) into (2.15) (more precisely, in Fourier 
transformation of (2.15), we define a well-known form [8, 22] current approximation which 
is linear with respect to a gauge field 

fiD.^k) = U>^%k)Al{k), (3.2) 

where 

p^(p-(fcgp) - ikp)d;)Kq 



U^"'{k) ^g^ J d^p 



pk + ipoe 

is the high temperature polarization tensor, and Meg = |(/° + /^) + Ncfg. 

Further we rewrite Eq. (2.16) in the momentum space. Taking into account (3.2), we 
obtain 

[A;V^ - (1 + Ok'-k'^ - n'^''{k)]Al\k) = jl'rik) - {jJI'm)+ 
+f''j s\;iTM{Al%k,)Al\k,) - {Al%kMl\k,)mk -k,- k,)dk,dk,+ (3.3) 
^f^cfff^ej Ea-t.3(^^^(^i)^I'(^2)A^1fc3) - (A^AkMl'ikMl'iks))) 

5{k — ki — k2 — ks) dkidk2dk3, 

where 

sltt = -^g{k^g''' + kW' - ki^g'^') , ^tt^ = aW^- (3.4) 

Let us multiply Eq. (3.3) by the complex conjugate amplitude A^*"(/c') and average 

it 

[k'g'"' - {l+r')k^k^ - ^'''{k)]{Al*\k')Al\k)) = {Al*'^{k')jlt{k))+ 

+f'"''l <K«*"(A:')^^^(A:i)>lI''(A:2))<5(A: - k, - k,)dk,dk,+ (3.5) 

^fcjjf,e j Eif^f^^,^{Al*'^{k')Al%^^^ -k,-k2- ks)dkAdks. 

We introduce the correlation function of the random oscillations 

Cik',k) = {Al*^{k')Al\k)). (3.6) 



In conditions of the stationary and homogeneous of QGP, when the correlation function 
(3.6) in the coordinate representation depends on the difference of coordinates and time 
AX — X' — X only, we have 

C{k',k) = C{k')S{k'-k). (3.7) 



By the effects of the nonhnear interaction of waves and particles, the state of QGP 
becomes weakly inhomogeneous and weakly nonstationary. The medium nongomogeneity 
and nonstationary lead to a delta-function broadering, and I^l depends on both argu- 
ments. 

Let us introduce I^ti^', k) = I^ti^, Ak), Ak ^ k' - k with \ Ak/k \<^ 1 and insert 
the correlation function in the Wigner form 

= J C{k,Ak)e-'^'^dAk, 

slowly depending on x. In Eq. (3.5) we change k ^ k' , a ^ b, complex conjugate 
and subtract obtained equation from Eq. (3.5), beforehand expanding of the polarization 
tenzor into Hermitian and anti-Hermitian parts 

W^ik) = U^'^^ik) + U^'"'{k) , U^'"'{k) = U*"'"'{k) , n^^^(A;) = -U*^'"'{k). 



We assume that anti-Hermitian part of 11"^ is small in comparison with 11^ and it is 
a value of the same smallness order, as the nonlinear terms in the r.h.s.. Therefore it 
can be suggested that 11^'"^ (/c) ~ 11^'"^ (/c'), and the term with can be rearranged to 
the r.h.s.. The remaining terms in the l.h.s. we expanded in a series in powers oi Ak to 
first smallness order. This corresponds to gradient expansion procedure usually used in 
derivation of kinetic equations [8, 12, 14]. Multiplying obtained equation by e~'^^^^ and 
integrating over Ak with regard to 



J Akx I^tik, Ak) Q-'^^^dAk 



dIfAKx) 



dx'^ 

we obtain finally 

— [k^g^'' - (1 + r')k^k^ - n^'^'^(^)]-^ = 2zn^'^^/;^- 



-z j dk'{{Al*'^{k')jlt{k)) - {Al\k)f^r{k'))}- 

-iif^ j dk'dk.dk, 45f£(^r"(^')^r(^i)>ir(^2))5(fc -k,- k,)- 

-r' j dk'dk.dk, ^,t2,X'(^r(^)^^*'^(^i)^r'(^2))5(^' -k,-k,)}- (3.8) 

.,{fcfff,ej 4f^';,\^,^{Al*^{k'^^^^^^ -k,-k,- k,)dk,dk,dk,- 



where jltik) ^ f^''^''^{k) + 

We make several remarks relative to obtained Eq. (3.8). The term with 11"^ introduc- 
ing in the r.h.s. of Eq. (3.8) corresponds to the linear Landau damping. However, as was 
shown by Heinz and Siemens [18], linear Landau damping for waves in QGP is absent 
and hence this term vanishes. Futher, the terms with S*^^) can be omitted also. These 
terms will be enter into kinetic equation for plasmons (Sec. 5) with Im E*^'^'' that vanishes 
by reality of the function E^^^. 

4. CONSISTENCY WITH GAUGE SYMMETRY 

In this Section we shall discuss the consistency of approximation scheme which we use 
with requirements of the non-Abelian gauge symmetry. 

The initial dynamical equations (2.4) and Yang-Mills equation (2.2) (without the 
gauge-fixing condition) transform covariantly under local transformations 

A^{X) = h{X){A^{X) + '-d^)h\X), h{X) = exp(^r(X)r), 

with parameter 9'^{X). We also have transformations of the quark, antiquark and gluon 
distribution functions [8] 

f,,^ip,X) = hiX)fg,g{p,X)h^iX), Up,X) = H{X)f,{p,X)H\X), 

where H"-\X) = ^^[t''h{X)t^h\X)]. 

As known (see, e.g. [16]), after the splitting of (2.5), (2.7) the resulting equations left 
two symmetries, the background gauge symmetry, 

= KX){A^^{X) + '-d,)h\Xl AliX) = h{X)Al{X)h\X), (4.1) 

and the fluctuation gauge symmetry, 

= 0, AliX) = h{X){A^^{X) + AliX) + '-d,)h'^{X). (4.2) 

The condition which we impose on regular part of a gauge field A^ and requirement 
that the statistical average of the fluctuation vanishes {A^) ~ 0, break both of types 
symmetry (4.1) and (4.2). Thus in the case of a gauge transformation (4.1) we obtain 
A^ 7^ 0, and in the case of (4.2) we come to noninvariance of the constraint {A^) = 0. 
Moreover, introduced correlation function (3.6) also have explicitly a gauge noncovariant 
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character. This leads to the fact that calculations in the preceding Sections are gauge 
noncovariant, and therefore their value is doubted. 

Nevertheless, there is the special case, when preceding (and following) conclusions are 
justified. This is the case of a colourless fluctuation, for which I^l{k,x) = 5°''^I^y{k,x). 
We can obtain a gauge invariant equation for I^uik, x) only in this restriction, in spite 
of the fact that the intermediate calculations spoil non-Abelian gauge symmetry of the 
initial equations (2. 2)- (2. 4). 

In principle, we shall be able to maintain explicit background gauge symmetry (4.1) at 
each step of our calculations, as it has been done, for example, by Blaizot and lancu [14] for 
derivation of the Boltzmann equation describing the relaxation of ultrasoft (A ~ 1/g'^T) 
colour excitations. First of all we assume that Ajj y^O. Then as the gauge-fixing condition 
for the random field Aj^, we choose the background field gauge 

V^{X)A''^{X) = 0, V^{X) ^d,- igA^iX), (4.3) 

which is manifestly covariant with respect to the gauge transformations of the background 
gauge field A^{X). We lastly define a gauge covariant Wigner function as in Refs. [8, 14] 

^) = / ^) s = X,-X2,x=^{X, + X,), 

where 

/J (., x) ^ C/-'(x, X + P'^x + |, ^ - |) U''\x - |, x), 

instead of the usual Wigner function /^^(/c, x), whose "poor" transformation properties 
follow from initial definition lf^{Xi,X2) = {Al\Xi)Al\X2)) . The function U{x,y) is 
the non-Abelian parallel transporter 

U{x,y)=Ve^v{-ig j^dz^'A^iz)}. 

The path 7 is the straight line joining x and y. 

The derivation of the kinetic equation for plasmons in this approach becomes quite 
cumbersome and non-trivial. For example, in the l.h.s. of equations for random parts 
of distributions (2.11)-(2.13), the covariant derivative will be used instead of the 
ordinary one d^. Besides, we cannot suppose that regular parts of distributions functions 
are specified and equal to equilibrium Fermi-Dirac and Bose- Einstein distributions (2.17). 
It is necessary also take into account their change using kinetic equations (2.8) with the 
collision terms in the r.h.s. of (2.8), which describe the reaction of the soft fluctuation on 
the background distributions. The correlators in the r.h.s. of Eq. (2.8) can be expressed 
in terms of the function /^^ and the distributions of hard particles f^ip, X), s — q,q,g , 
only. 

However, if we restrict our consideration to study of colourless excitations and replace 
the distribution functions of hard particles by their equilibrium values (2.17), then it 
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leads to effective vanishing of terms explicitly contained the mean field Aj^. This follows 
e.g., from the analysis of derivation of Boltzmann equation by Blaizot and lancu [14]. 
Therefore the most simple way of derivation of kinetic equations for soft colourless QGP 
excitations is assume = and use prime gauge noncovariant correlator (3.6). In this 
case the background field gauge (4.3) is reduced to a covariant one. The kinetic equation 
obtained by this means is gauge invariant if all contributions at the leading order in g 
to the probability of nonlinear wave scattering by hard thermal particles are taken into 
account (see also discussion in Conclusion). 

In fact, the requirement of nontriviality of a colour structure of Wigner function 
I^^{k,x) leads to necessity of existence of nonvanishing mean field A^ and/or the ex- 
ternal colour current and conversely, the existence of mean field in QGP results in colour 
structure I^l{k,x), which is different from the identity. 



5. THE SECOND AND THIRD APPROXIMATIONS OF THE COLOUR 
CURRENT 



Now wc conscrncd with computation of the nonlinear corrections to the current in the 
r.h.s. of basic Eq. (3.8). 

At first we define /^^^'. We carry out the Pourier transformation of Eq. (2.12) and 
substituting the obtained /^^^^ from (3.1) into derived expression, we find 

pK -\- ipoe J pK2 -T ipo€ 

(5.1) 



5{k — ki — k-i) dkidk2- 

Prom here on the suffix "T" for a gauge field is omitted. The expression for fj^'^\k,p) 
is obtained from (5.1) by choosing upper sign and replacements fg^^"^ ~^ 

Substituting obtained expressions /J^^-*, s — q,q,g into (2.15) (for n = 2), we find required 
current correction 

^•T(2)a.(^) = -^g^r^'■ f (^M^ (^^(^^1) ^^(^2) - {Al{k,) Alik2))) (5.2) 

J pk + ipoe pk2 + ^poe 

8{k — ki — k2) dkidk2+ 

4-- J ^^^-^^ ^) (-'(^.).^(.)-<.t(..)<(.,)) 

5{k — ki — k-i) dkidk2- 
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The contribution of gluons to the expression with symmetric structure constant d"'^'^ here 
drops out. This is connected with the fact that in calculation of trace of anti-commutators 
we have: Spr{^^^^} = ^0?"'"= - for quarks and antiquarks, and Tr^"{^^^^} = - for 
gluons. The symmetry of contributions can be restored if we note that besides usual gluon 
current jg{x) = gt"- j d^pp^^TY{T"'fg{x,p)), the dynamical equation for gluons admits a 
covariant conserving quantity 

Cgf^l d^pTV(P%(^,p)), (5.3) 

where {V"-)^'^ = d"''"^ and is a certain arbitrary constant. The covariant continuity of 
(5.3) is evident from the identity: [P^T^] = iP^^V". On addition of (5.3) to (2.3) we 
have contributions to the nonlinear current corrections only. Adding (5.3) to the second 
current iteration (2.15) and taking into account the equality 

we derive more general expression for instead of (5.2) 

^■T(2)a,(^) = j 5$r,J(4(^l)^^(^2) - {AlikMlihWik -h- k,)dhdk„ (5.4) 

Where ^^^^^ = /'^'^^S^' + d^'^SifX^ 



' ' J pk + ipoe pk2 + ipoe 

^kMM 2 J ^ pk + ipoe dp- \pk2 + ipoe dpP ) ' ^ ' 

1 

2' 



The tensor structure of S^^lp^^^ exactly coincides with appropriate expression ob- 
tained in calculation of j'^*^^-*'^ in Abelian plasma [21], and hence piece of current with d"''"^ 
has a meaning of Abelian part of the colour current xhe term with S\!J^\^I^^ is 

purely non-Abelian, i.e. it has no Abelian counterpart. 

Let us estimate orders of 5'^-'^^^ and S'^^^^\ Following usual terminology [1], we call an 
energy or a momentum "soft" when it is of order gT, and "hard" when it is of order T. 
We will be considered, as in Ref. [12], that collective excitations carrying soft momenta, 
i.e. k ~ gT, and plasma particles have the typical hard energies: p T. In a coordinate 
representation the first of conditions denotes that oscillation amplitude A^^ and distribu- 
tion functions of hard particles fs{X,p),s = q,q,g change on the scale X ~ i/gT. On 
this basis, we have the following estimate for 5'^^''^^ 

slfX - g'T. (5.7) 
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Here, we considered that by virtue of the definitions (2.17) Meq ~ 1/T^. 

In expression (5.6) the integral of energy with gluon distribution function is infrared 
divergent. In a similar manner [17], we regulate it by introducing an electric mass cutoff 
of order gT, and only take the leading term in g. Then it can be found that in (5.6), the 
part related to the gluon distribution function is of order g^T{g In g) and the other part 
related to the quark and antiquark distribution functions is of order g^T, i.e. 

s\:ZC-9'ng\^g)+g'T. 



Hence, 5'*^^^-' which is purely non-Abelian, is of lower order in the coupling constant 
than S'^^^^\ which has an Abelian counterpart. This fact was first seen in Ref. [17]. 

The expression for a colour current in third order in field is defined by means of 
reasoning similar previous ones. Performing the Fourier transformation of equation (2.13), 
taking into account (3.1), (5.1) and equalities 

Tr ({T", T^jfT"^, T"}) = NJ'^V^ + 46''^6"^ + 25"'^5"'' + 25"^^'"^, (5.8) 

Tr({7'",r''}{r'^,r^}) = o, 

we find the required form of j^(^)"'* 

-{Al{h)A'i{h)Al{h)))5{k -h-h- h) dhdk2dh+ (5.9) 

+da,cfc,el i?^;;-;^,^^,^(A^(A;3)A^(A;0A^(^2) - {Al{ks) A'i{k^) AUk^))) 

S{k — ki — k2 — k^) dk\dk2dkz. 

Here, 

^kMMM ~ / / ^kMMM ^ J " ^kMMM " " ^kMMM ^ J ^kMMM'^ 

kMMM y J ^ pk + ip^e p{ki + k2) + ipoe p/ca + ipoe ^ ' ' 

kMMM J ^ pk + ipoe dp"- ^p{ki + ^2) + ipoe dp"' ^pk2 + ipoe dpp 

(5.11) 

yiVI)^luXo■ _ (IV)^uXa 
k,ki,k2,kz 2 ^^^\ik2,kz' 

The expression for S*^^^^) is obtained from (5.11) by exception of quark and antiquark 
contributions. The availability of the term with YP'^^'> is refiection of more complicated 
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colour structure of the gluon kinetic equation in comparison with quark and antiquark 
equations, that manifests here, in appearance of additional terms in (5.8) as compared 
with 

Sp({i", t''}{t'^, t'}) = ^rt"'^" + -^S'^'^S'^. 

The terms with E^^^^), E^^) and R are defined as the interference of Abelian and non- 
Abelian contributions. For colourless fluctuations of QGP, which we study in a given 
paper, the correlation function (3.6) is proportional to the identity. This leads to the 
absence of the interference of Abelian and non- Abelian contributions. For this reason 
their explicit form is not given here. 

At the end of this Section we estimate the order of S*^'^^-' and S^^^^. It follows from 
the expression (5.1) that 

SgSa - (5-12) 

Cutting off, as in the previous Section, integration limit for the gluon distribution function, 
we flnd 

E(^^)~E(^^)~y3^/, E(^^^)~/. 

By this means, purely non- Abelian contribution of E^^^^ is of lower order in the coupling 
constant than Abehan - E^-^^), E(^^) and E^^^^). 



6. THE KINETIC EQUATION FOR LONGITUDINAL WAVES 

Now we turn to initial equation for waves (3.8). We substitute obtained nonlinear 
corrections of induced current by field (5.4) and (5.9) into this equation. 

Because of a nonlinear wave interaction the phase correlation effects take plays. By 
virtue of their smallness, fourth-order correlators can be approximately divided into prod- 
uct of the correlation functions {A* {k') A(k)) . For third-order correlation functions this 
decomposition vanishes, and it should be considered a weak correlation of phases fields. 
For this purpose we use the nonlinear equation of a field (3.3), taking into account only 
the terms of the second order in A, in the r.h.s. of (3.3) 

[ytV - {l+C')k''k'' - ni"'{k)]Al{k) = 

= r'^/ Sj::,l,,^{Al{kMl{k2) - {Al{k-,)Al{k^)))5{k -k,- k^)dk,dk^. (6.1) 
The approximate solution of this equation is in the form 

Alik) = Af'^{k)-V,^{k)r"^ j S'^},l,^{Af\kM^Hk2)- 

-{Af\ki)Af\k2)))6{k -ki- k2)dk,dk2, (6.2) 
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where S^f^^^^ = + Si^MM^ ^f^"!^) is a solution of homogeneous Eq. (6.1) 

corresponding free fields, and 

T>,,Ak) = -[k^Qi.. - (1 + r')ki,K - n^.(^)]-' (6.3) 

represents the medium modified (retarded) gluon propagator. 

Now we substitute (6.2) into third-order correlation functions entering to Eq. (3.8). 
Because A^^^ represents amphtudes of entirely uncorrelated waves, the correlation func- 
tion with three A^^^ drops out. In this case every term in {A*^(k')Al(ki)Af(k2)) and 
{A^^{k)Al'^{ki)A*^{k2)) should be defined more exactly. In the correlation functions of 
four amplitudes, within the accepted accuracy, it can not make distinctions between the 
fields A and A^o). 

Finally Eq. (3.8) becomes 

Q gjab 
-[A;V" - (1 + r^)k''k'' - n^'^^(A;)]- - 



dk\ dx^ 

^-ij dk'dk,dk,dk,{tfff'^5{k-k,-k2-k,)t';:^^^^^^^^^^ 

-{Al\k')Al{k,)){Ai{kMl{k2)))- 
-r^ff'^6ik' -k,-k,- A:3)S*r,^:,,,,3((4(A:)Ar(A;3)A*'^(fci)Ar(fc2))- 

-{A%k)A:\k,)){Al\k,)Al\km+ (6.4) 

^ifcdfaefj j dk,dk2dk[dk',{v;^{k')-v^,{k))si^;:^^^^^^ 

Ank[)A*J{k',)) - {Al{kMiik2)){ATik[)A*J{k',)mk -k,- k2)6{k' -k[- k',). 
Here, we keep the terms of leading order in g only and set 

It follows from the definition (6.3) that the propagator is of the order ~ 
Taking into account (3.4), (5.7) and (5.12), we see that all terms in the r.h.s. of (6.4) are 
of the same order. This explains the fact that in the expansion of the current (2.14) the 
following term - j^*^^-* should be retained in addition to the first nonlinear correction j^*^^^: 
it leads to the effects of the same order of magnitude. 

Let us make the correlation decoupling of the fourth-order correlators in the r.h.s. of 
Eq. (6.4) in the terms of the pair ones by the rule 

{A{k,)A{k2)A{k^)A{ki)) = {A{k,)A{k2)){A{k^)A{ki)) + {A{k,)A{k^)){A{k2)A{ki))+ 

+{A{k^)A{k^)){A{k2)A{k^)). 
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Taking into account that the spectral densities in the r.h.s. of Eq. (6.4) can be 
considered as stationary and homogeneous those, i.e. having the form (3.7), and setting 
= (5»&7^^, we find instead of Eq. (6.4) 

r) f)T 
-[fcV" - (1 + i~^)k^k'' - n^'''^(A;)]^ - 



dkx dx^ 
= I dk, lmit>^,%^_,^ - t'i^l\_,;)I,x{k)I.,{k,)+ (6.6) 
+NM{V^{k)) I dkA {S^,%1,^ - Si%^,;){SlfX - Sl%>:,,)I,>.{k,)hAk,)5{k -k,- k,). 



As it is known [23, 24], in global equilibrium QGP the oscillations of three typies can 
be extended: the longitudinal, transverse and nonphysical (4-D longitudinal) ones. In 
this connection we define the Wigner function I^i,{k, x) — I^v in the form of expansion 

V = P^^A + Q^^A + iD,Jk , 4''^'"^ = I^^'^'^^Kk, x). (6.7) 
The Lorentz matrices in (6.7) are members of the basis [24, 25] 

Pnv{k) — Qfj^i, — Dn^(k) — Qfj,u(k), Qij,iy(k) — -^^4y7Tr^) Cfj_i,{k) — — ^'^^^ ^ ) — 

D^u = k^ky/k , u^{Jt) = k k^iykxi). 
The effective gluon propagator (6.3) can be written in more conveniet form 

V^,{k) = -P^,{k)^\k) - Q^,{k)^\k) + iD^,{k)^\k), (6.8) 

where A*''(A;) = 1/(A;2 - n*''(A;)), n* = P^^.Jl^ ^W^Q^y-^^{k) = 1 / {{u + ief - k'^) . 
The shift ie is introduced in A^{k) to provide the right analytical properties. At finite 
temperature, the velocity of plasma introduces a preferred direction in space-time which 
breaks manifest Lorentz invariance. Let us assume that we are in the rest frame of the 
heat bath, so that = (1, 0, 0, 0). 

Further derivation of kinetic equation for longitudinal oscillations involves the same 
type of manipulations as in the theory of electromagnetic plasma, so we can afford to be 
sketchy. 

Now we omit nonlinear terms and ant i- Hermit ian part of the polarization tensor in 
Eq. (3.5). Further substituting the function 5°'^Q^y{k)ll.5{k' — k) instead of I^t{k', k), we 
lead to the equation 

Re{e\k))li^O. 

Here, we use relation 

A-^'(A;) = A;V(A;), (6.9) 
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where 

l + x 



l-x 



m9{l - \x\)] 



is the longitudinal colour electric permeabihty and to'^i = j^g'^T'^{Nf + 2Nc) is a plasma 
frequency. The solution of this equation has the structure 

li = Il6{u - ool) + I'_^5{u + ul,) , 4 > 0, (6.10) 

where is a certain function of a wave vector k and a;^ = a;'(k) is a frequency of the 
longitudinal eigenwaves in QGP. 

The equation describing the variation of spectral density of longitudinal oscillations is 
obtained from Eq. (6.6) by replacement: I^i, — > Q^v{k)l]^, where l\. is defined by (6.10). 
5-functions in (6.10) enable us to remove integration over frequency and thus we have 
instead of Eq. (6.6) 

/ ^ 2^^^, I rfk,/i.(im[(Ej£;.,_„ - EKl-*.)+ 

+ (6.11) 

roc /■ ^ ^ 



where 

GkMM = Mi^pa{k)KS'k%:M - S^mm){SI%\% - SljX)Q,>^ik,)Q.<^{k,)5{k -k,- k,). 

Let us consider the terms in the r.h.s. of Eq. (6.11). The integral with the function 
jfcj^jfcj is different from zero if the conservation laws are obeyed 

k = ki + ks, 

These conservation laws describe a decay of one longitudinal wave into two longitudinal 
waves. However for a dispersion law of the longitudinal oscillations in QGP, these reso- 
nance equations have no solutions, i.e. this nonlinear process is forbidded. Therefore the 
integral with Gk,kiM vanishes. Remaining integrals with G-functions differ from Gk,kiM 
in that some of the interacting waves are not radiated but absorped. They also vanish. 

The expression 

(^fe,jk,-iti,fei ~ ^fc,-fci,fe,iti) Q/iiA(^)Qi/<T(^l) |a;=w|^,a;i=a;J^^) (6-12) 

contains the factors 

1/ {pk + ipoe) , 1/ {pki + ipoe) , 1/ {p{k - ki) + ipoe). 
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Imaginary parts of first two factors should be setting equal to zero, because they are 
connected with linear Landau damping of longitudinal waves (which is absent in QGP), 
and therefore, the imaginary part of the expression (6.12) will be defined as 

= -jj^5(4-<-v(k-kO). 

It follows that nonlinear term in the r.h.s. of (6.11) with the function (6.12) is different 
from zero if the conservation law is obeyed 

4-<-v(k-ki) = 0. (6.13) 

This conservation law describes the process of scattering of plasmon by the hard thermal 
particle in QGP. 

Let us consider in more detail the term in (6.12) (sec definition (6.5) with propagator 
T^paik — ki). By expansion (6.8) this propagator represents the nonlinear interaction of 
longitudinal waves with longitudinal ones through three types of intermediate oscillations: 
the transverse, longitudinal and nonphysical oscillations depending on a gauge parameter. 
The term with 

in general, describes two fundamentally different nonlinear processes: 

1. if (cUk — uj\^^, k — ki) is a solution of the dispersion equation A*(/c — /ci) — 0, then 
this term describes the fusion process of two longitudinal oscillations in transverse 
eigenwave; 

2. otherwise, it defines the process of nonlinear scattering of longitudinal waves in 
longitudinal those through the transverse virtual oscillation (for a virtual wave in 
distinction to the eigenwave, a frequency uj and a wave vector k are not connected 
with each other by the dispersion dependence 7^ tc;(k)). 

The equahty A'(/c — ki)\^^^i ^,=1^1 — does not hold for longitudinal oscillations, as we 
see above and therefore, the term 

(Q,aik-k,)A\k-ki)) , 

defines only the process of scattering of longitudinal waves in longitudinal those through 
the longitudinal virtual oscillation. 

The contribution of nonphysical intermediate oscillations 

UD,^{k-k,)A%k-k^)) ^ 



Im 



1 

vik -kA + ivne 
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will be considered in Sec. 12. 



The remaining terms in the Eq. (6.11) with S are distinguished from above considered 
ones by a sign of ki, and describe the processes of simultaneous radiation or absorption 

by particles of two plasmons. 



Summing the preceding and going to the function 

,dRee\k)\ 



Wi 



dijj 



-'k > 



having the physical meaning of spectral density of longitudinal oscillations energy, we find 
from (6.11) the required kinetic equation for longitudinal waves in QGP 



where 



Vl 



dujl 
dk 



Wl 



dRee^{k). /.dRee\k) 



(6.14) 



dk 



duj 



is the group velocity of longitudinal oscillations and 

■WL^\ 1 (dRee\k)^-i(dRee\ki)^-i 



k 



y(k)=27Ve/ rfki(-J 



/crtefc \^v)\-^ /usx^t \rii)\- 



I™[(^fc,fc,-fci,fci ^fc,-fci,fc,A;i) + (^fc,fe,fci,-fci ^kM,k,-k])\Q h\{^)Qi"t{^i) 



(6.15) 



presents the damping rate caused by nonlinear effects and being the linear functional of 
spectral density of energy. 

One can write (6.14) in more convenient form if the spectral density of number of 
longitudinal oscillations is represented as 

It fulfils the role of distribution function of plasmons number. Then instead of (6.14) we 
have 

dNl _ dNl dNi 



dt 



dt 



ax 



-7 {K} K ■ 



(6.16) 



7. HTL- AMPLITUDES. WARD IDENTITIES 



Befor proceeding to an queshion on a gauge dependence of obtained nonlinear Landau 
damping rate (6.15), we rewrite derived expression in the terms of HTL-amplitudes [1, 
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12]. This makes possible to extend a procedure of a gauge invariance proof of the damping 
rate for QGP collective excitations in quantum theory [1] to our case. 

We present the integration measure d'^p as dp°|ppd|p|dr2, where dQ is the angular 
measure. Using the definition of equilibrium distributions (2.17) (for = 0) and taking 
into account 

—oo —oo 

we perform the integral over dpo and the radial integral over d|p| in the expressions for 
5'(-f^)_function (5.5) and E^^^^-function (5.10). Futher we rewrite the expression (6.15) as 
follows 



(k) = -2/iv./dk,(-i)[j5j,(^-U) ]_,^tor(k,k,), 

(7.1) 

where 

T(k,ki)^ (7.2) 
= {6r''''^''{k, ki, -k, -ki) - *r^"'P{k, -ki, -k + ki)Vp^{k - ki) T"^"(A; - ki, -k, ki) - 
- *r^'"'{k, k,, -k - k^)Vp^{k + k^) T"^'^(A; + k,, -k, -k,)}Q^x{k)Q,a{h) 

and 

dr {k,k„k„ks)-3u;^,j + + 

' ^ '\{v^^il,.)) (7.3) 



vk^~ie^ v{k + k-^) + ie^vki — ie vk2 — ie 
present HTL-corrections to the bare four-gluon vertex [1, 12]. 

*r^"'P{k, ki, k2) = V^^^ik, ki, k2) + ST^^'Pik, ki, k2) (7.4) 

is the effective three-gluon vertex [1, 12], that is a sum of the bare three-gluon vertex 

V^^'^'ik, k,, k2) = g^'^ik - kiY + g^'^iki - k2Y + 9^''{k2 - k)" (7.5) 

and corresponding HTL-correction 

4„*2) = 3a,J / ^ {-^ - -^). (7.6) 

' ' P^J A-K vk + ie^vk2-ie vk^-ie^ ^ ' 

The polarization tensor in these notations takes the form 



dn yf^v" 
An vk + ie 
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For writing the expression (7.1), for example, in the temporal gauge Aq = 0, it is 
sufficiently to replace the projection operators Q^\{k) and Qua{ki) by 



(similarly for Qvaiki)), and the propagator (6.8) by 

-^D,^(k)A'(k) - ^,-^D,^{k), (7.7) 

where is a gauge parameter in the temporal gauge. 

To establish the gauge invariance of nonlinear Landau damping rate 7'(k) there is a 
need to show that expression ImT'(k, ki), where function T'(k, ki) defined by (7.2) in a 
covariant gauge equals to similar expression in the temporal gauge. 

We prove a gauge invariance for more general expression that in a covariant gauge has 
the form 

*r{k,-k2,k,,-k3)= (7.8) 

= { T^^^'^(A;, -k2, ki, -ks)-V,a{-ki+k2) *r'''''{k, -ks, ki-k2) T"^'^(-A;i+A;2, k^, -k2)- 

-Vpai-ki+ks) *r>"'P{k, -k2, ki-ka) T"^'^(-A;i+A;3, ki, -k3)}u^{k)ux{ki)uMMk2)\on-sheii. 

Here, k + ki = k2 + k^. The effective vertex *Yi^'^^^ are formed by adding the hard thermal 
loop (7.3) to the bare four-gluon vertex 

The association of the expression (7.2) with (7.8) is given by 

Imf(k,ki) = \ Im*f(fc, -k2, fci, -/i:3)|fci=-fc, fc2=-fci, fc3=fci- 
w^lkjw^lki) 

As will be shown in a separate publication [26], the expression (7.8) is associated with 
probability of plasmon-plasmon scattering. 

Similar expression (7.8) in the temporal gauge is obtained with replacements Un{k) — > 
Un{k) and the propagator (6.8) by (7.7). 

The gauge invariance proof is based on using the identities, analogous the effective 
Ward identities in hot gauge theories [1]. It can be shown that the following equahties 
hold 

k^*^f''^^-^k,k^,k2,ks)= *T''^''{ki,k2,k + ks)- *T^^''{k + k^MM), 
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k,^*Y'^'^^-{k,h,k2,h)= *T'^^''{k + h,k2,h)- *r''^''{k,ki + k2,k3), (7.9) 
(similar contractions with k2x, k^^), 

k^*T''''''{kMM) ^'D-^'^'i-ki) - V-^^p{-k2), (7.10) 

(similar contractions with kiy,k2p). 

Here V-^t"'{k) = P^"'{k)A-^\k)+Ql"'{k)A-^^{k) is the inverse propagator for which 
the following useful relation holds 

T^paT^-"'\k)=S'^-^. (7.11) 



8. THE GAUGE INVARIANCE 



Our proof of gauge invariance is reduced to contraction of the projectors Q and Q 
with resummed three- and four-gluon vertices and the use of the Ward identities (7.9), 
(7.10). At the beginning we consider the term with a gauge parameter for the propagator 
in a covariant gauge. 

By using the Ward identities (7.10), we have 

^D,^ik, - k2)A%k, - k2) *T^''^{k^, -ks, k, - k2) T"^'^(-A;i + ^2, k^, -k2) = 

= e(A°(A;i - k2)f{V-^^''{k) - V-^''''{ks)){V'"'{k2) - V-'^^i-k,)). (8.1) 

Further this expression is contracted with Ui^{k)ui,{k3). It is easily shown that it vanish 
whether because 'D~^^^{k) is transverse, or by the definition of the mass-shall condition, 
i.e. 

k,V-^^\k) = 0, I?-^°^(A;)L=,z = 0. (8.2) 

Similar statement holds in the temporal gauge also. 

The gauge- dependent parts in the above calculation drops out 7'(k), since they are 
multiplied by the mass-shell factor. These factors are proportional to (u — ujIJ. However, 
in a quantum case Baier, Kunstatter, and Schiff [27] observed that naive calculation in 
covariant gauge appears to violate this consideration. Mass-shell factor is multiplied by 
the integral involving a power infrared divergence which is cutoff exactly on the scale 
(a; — 0;^) ~ g^T. By this means the gauge-dependent part yields a finite contribution to 
the gluon damping rate. This problem is considered for damping rate of Fermi-excitations 
in QGP, also [28, 5, 29]. 

It can be shown that in our case, the integral preceding the mass-shell conditions, 
diverges for lower limit also and thus the similar problem is arised: does (8.1) yield a 
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finite, the gauge-dependent contribution to the nonhnear Landau damping rate? In Sec. 
12 we provide the answer to the queshion. 

Now we consider the remaining terms in (7.8). We calculate the contraction with 
effective four-gluon vertex *r4. Slightly cumbersome, but not complicated computations 
by using the effective Ward identities (7.9), (7.10) and relations (7.11), (8.2) lead to the 
following expression 

*r'^^''{k, -k2, h, -h)u^{k)Ux{ki)u^{h)u„{k2)\on-shell = 

= k'^klklkl T°°°°(A;, -k2, ki, -ks) + E{k, -k2, ki, -ks), (8.3) 

where 

Eik, -k2, ki, -ks) = {{klkliukl + u^k^] T°°°(A; - ks, ki, -ki)- 
-uusklklksu "V^ik - ks, ki, -/ca)) + {k ^ ks) + {ki ^ k) + {ki ^k,k2^ ks)}+ 
+{{{ujkl + usk''){uj2kl+ujikl)V-^''''{-ki + k2)-ujiUJ2iujkl + ujsk^^^^ 
-ujujsiuj2kl + uikDks^V-^^^i-ki + k2) + ujuJiUJ2UJski\ksv'D'^''\-ki + /cs)) + (/ca ^ h)}- 

Further, we calculate the contractions with the terms containing in (7.8). Here, 
we are led to the expression 

{Dpa{-ki + k2) *r'"'P{k, -ks, ki - k2) T"^'^(-A;i + k2, k^, -k2) + {k2 ^ ks)}u^{k)ux{k{) 

Uu{ks)u„{k2)\on-shell = (8.4) 

= k^klkjkliVp^i-k, + k2) *T^'P{k, -ks, ki - k2) T"°0(-A;i + k2, k,, -^2) + (^2 ^ ks)} 

+E{k, -k2, ki, -ks). 

Note that the function S in the r.h.s. of (8.4) is identical to that in (8.3). Subtracting 
(8.4) from (8.3), we arrive at desired expression 

*rik, -k2, ki, -ks) = k^klklkU T°°°°(A;, -k2, ki, -ks)-Vp^{k-ks) T°'^''(A;, -A;3, ki-k2) 
T"00(-A;i+A;2, fci, -k2)-V,^{k-k2) T0°^(A;, -k2, k,-ks) T"00(-A;i+A;3, k,, -ks)}. (8.5) 

Now we consider the structure of *r in the temporal gauge. For this purpose we 
replace -u^ by -u^ in (7.8) and the propagator in the covariant gauge by the propagator in 
the temporal gauge (7.7). The contraction with effective four-gluon vertex leads to 

*r'^''^''{k, -k2, ki, -ks)Ui_i{k)ux{ki)u^{ks)Ua{k2)\on-shell = 

= eklklkU T°°°°(A;, -k2, ki, -ks) + E{k, -k2, ki, -ks)}, (8.6) 

where 

E{k, -k2, ki, -ks) = {{luilo2{lo + Us) T°°°(A; - ks, k^, -k2)- 
-ujiuj2ks^ T''°°(A; - ks, ki, -k2)) + {k ^ ks) + {ki ^ k) + {ki ^k,k2^ ks)}+ 
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+ u;2){uj + cc;3)p-^°°(-fci + ^2) - {u + U3)kixV-''>\-ki + ^2)- 

-{ui + u2)h,v-^''''{-h + ^2) + hxh.v-^''\-h + ^2)) + (^2 ^ ^3)}- 

Contracting with the terms containing *r3, yields 
{Vpa{-ki + k2) *V^'"'{k, -k^, ki - k2) *r"^''{-ki + k2, ki, -k2) + {k2 ^ k^)}u^{k)ux{ki) 

Uu{k'i)u„{k2)\on-shell ^ (8.7) 

= eklklkl{Vpa{-ki + k2) T°°^(A;, -A;3, ki - ^2) T"°°(-A;i + k2, /ci, -A;2) + (^2 ^ ^3)+ 

+E;(/c, -A;2, /ci, -A;3)}. 

Subtracting (8.7) from (8.6), we are led to similar expession (8.5). Thus, we have shown 
that at least in a class of the covariant and temporal gauges, the nonlinear Landau damp- 
ing rate (7.1) (exactly, its piece, independent of a gauge parameter) is a gauge- invariant. 

9. THE PHYSICAL MECHANISM OF NONLINEAR SCATTERING OF 
WAVES 

In Sec. 7 the expression of the nonlinear Landau damping rate 7'(k) for Bose excita- 
tions in a quark-gluon plasma was obtained. Now we transform it to the form, in which 
each term in 7'(k) has direct physical relevance. The first transformation of this type has 
been proposed by Tsytovich for Abelian plasma [21]. In [20] it taken into account the 
contribution with longitudinal virtual wave only. In order to satisfy gauge invariance, we 
extend the transformation of similar type to the case with the transverse virtual wave. 
Since the nonlinear Landau damping rate is independent of the choice of gauge, we choose 
the temporal gauge for simplicity. 

As was mentioned in Sec. 6, the expression for 7'(k) contains the contributions of two 
different processes. The first is associated with absorbtion of a plasmon by QGP particles 
with frequency u and a wave vector k with its consequent radiation with frequency ui 
and a wave vector ki. It defined by the second term in the r.h.s. of (7.2). The frequency 
and wave vectors of a incident plasmon and a recoil plasmon satisfy the conservation 
law (6.13). The second process represents simultaneous radiation (or absorbtion) of two 
plasmons with frequency a;, uji and wave vectors k, ki satisfying the conservation law 

4 + <-v(k + ki) = 0, (9.1) 

and it defined by third term in (7.2). In contrast to previous scattering process, this 
process not conserves the plasmons number and its contribution are not important to the 
order of interest. The first term in (7.2), associated with HTL-correction to four-gluon 
vertex, contains both processes. Futher, we have taken into account the terms with (6.13) 
and we drop the terms which contain a (5-function of (9.1). 
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Now we consider the expression with 6T4. With regard to above-mentioned and by 
the definition (7.3), the contribution of a given term to 7'(k) can be represented as 

f /W^ \ dO 

3MgfJ (4-O<(k,k0 l^-^j —dk„ (9.2) 



where 



<(k, ki) 

To clear up the physical origin of contribution (9.2), it is convenient to compare it with 
corresponding contribution in the theory of electromagnetic plasma. In this case as shown 
in [21] this contribution describes the normal Thomson scattering of a wave by particles: 
a wave with the original frequency is set in oscillatory particle motion of plasma 
and oscillating particle radiates a wave with modified frequency o;^^. The corresponding 
function w^{]<., ki) presents the probability of Thomson scattering. As it was shown above 
in quark-gluon plasma for a soft long-wavelength excitations all Abelian contributions is 
at most glng times the non- Abelian ones and the basic scattering mechanism here, is 
essentially another (our discussion is schematic, the details of a computation displayed in 
Appendix A). 

To show this mechanism we use the classical pattern of QGP description [8], in which 
the particles states are characterized besides position and momentum, by the color vector 
(non- Abelian charge) Q = (Q"), a = 1, . . . , — 1 also. As was shown by Heinz [8] there 
is an close connection between the classical kinetic equations and semiclassical ones (2.4). 
Therefore in this case the use of classical notions is justified. 

Let the field acting on a colour particle in QGP represents a bundle of longitudinal 
plane waves 

A^^ix) = -J [(a.VA;^)Q,.(A:)Ar].=.j^e^'^---l*dk. (9.5) 
The particle motion in this wave field is described by the system of equations 

= SQ'F"^'^, (9.6) 

= -sf-^Alq'. (9.7) 

Here, r is a proper time of a particle. Eq. (9.7) is familiar Wong equation [30]. The system 
(9.6), (9.7) is solved by the approximation scheme method - the weak field expansion. A 
zeroth approximation describes uniform restlinear motion, and the next one - constrained 
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charge oscillations in the field (9.5). With a knowledge of the motion law of a charge, the 
radiation intensity by it longitudinal waves can be defined. In this case Eq. (9.6) defines 
the Abelian contribution to radiation, whereas (9.7) - non-Abelian one and interference of 
these two contributions equals zero. The scattering probability computed by this means, 
based on Eq. (9.7) is coincident with obtained above (9.3). 

In this manner the contribution of (9.2) to 7'(k) are caused by not the spatially 
oscillations of a colour particle, as it occurs in electromagnetic plasma, but it induced 
by a precession of a colour vector Q of a particle in field of a longitudinal wave (9.5) 
(Q«Q" = const by Eq. (9.7). 

Let us consider now more complicated terms in (7.1) associated with * Fa-functions. 
By the definition of three-gluon vertex the following equality is obeyed 

*r''''P{k, -h, -k^) = - T'"^'^(A;2, -K h), 

(hereafter, k2 = k — ki). Using this relation, the contribution of the term with longitudinal 
virtual wave to 7'(k) can be presented as 



-2/7V, 



1 



Futher we use the relation 



y k^kf kl ^ duo Ji^=ujl\ duji ^'^i='^Li 
Im {kll\\k2){*r^\k, -ki, -k2)k'k{k^^f) 



a)=a){^, 0)1=0) 



l^{kl/\\k2){*V^^k'k{k},f) = -\m{kl^\k2))-^\kl^\k2){*r^'k'k{k'2)\^+ 

+2(Im *r^^k'k{kl^)Re{kll^\k2){*r^^k'k{kl^)). (9.8) 

Here, we drop the momentum dependence on vertices. Taking into consideration the 
equality 

Iu.{kl^\k2))-'^2>^ul'^j'^5{vk2). 

one can write contribution of the first term in the r.h.s. of (9.8) to the nonlinear damping 
rate in the form 



r /W'- \ HO 



where 

w|(k,ki) = 



k ^ — ^ ■ '^i='^k 
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MHK^i) = Tj^^^(^^^^'''ik,-h,-k,)k%kl^ . (9.9) 

|k| |ki| ki V ^^i^i /.=.j^,.,=4 



HTL-correction ST^ enters in the amplitude (9.9) only. Here the contribution of bare three- 
gluon vertex is dropped out, since by the definition (7.5), P^'(A;, -ki, -k2)k'k{k^^ = 0. 

The quantity Ai " can be interpreted as the scattering amplitude of a longitudinal wave 
by dressing "cloud" of a particle. However, in contrast to the Abelian plasma here, the 
scattering is produced not by the oscillation of a screening cloud of color charge as a result 
of interaction with incident scattering wave, but as a consequence the fact that it induced 
by precession of color vectors of particles forming this cloud in the incident wave field 
ul-. This process of transition scattering is a completely collective effect. For calculation 
of its probability tt;J|(k, ki) it is necessary to solve the kinetic equation describing a color 
charges motion in a screening cloud in the field that is equal to the sum of fields of incident 
wave (9.5) and a "central" charge producing a screening cloud. 

Now we consider the contribution with transverse virtual wave. Using the association 
of projectors 

^2 ^ 

Ppaik'l) = 9 pa ~ UpUa ~ ~r2Qpa{k2), 

we rewrite it in the form 

2a'N fdk (^^A ^ ( dRes\k) yi . dRes^(k,) yi 
-2g N^J dki l^-^J k%!l^^J.=.<l Ou;, i.,=4, 

Im (A*(A;2) *r^\k, -k^, -k2)k% *r'^'\k, -ki, -k2)k''k(- (9.10) 



~^\k2){*r'\k, -k,,-k2)k%kif] 



To simplify the expression, we expand *V^^^k'^k[ in two mutually ortogonal vectors: \^2 
and [nk2], where n = [kki] 



hi [rib-, 1^ 

r'^^k% = ( *r^'k%k'2)T\ + ( *r^'k%[n\i2Y) ^ 



Substituting the last expresion into (9.10), instead of the expression in square brackets, 
we define 



^ ^ I^{i^\k2){*r^\k,-ki,-k2)k%[n\i2]'f) 



'^=4 '^1=4/ 



We transform the imaginary part similar to (9.8) 

Im(A*(A;2)(*P^'A;^A;i[nk2]0') = -Im(A-^*(A;2))| A*(A;2) *P^Wi[nk2]f + 

+2(Im*rJ'A;^A;{[nk2]0Re(A*(A;2) T^^k'k{[n\i2]^). (9.11) 



28 



If wc take into account Im A ^*(A;2) = — (l/2)ImA ^'(^2), then the contribution of the 
first term in the r.h.s. of (9.11) to 7'(k) can be presented in the form 



SMgfJ (4-O«^v^(k,k0(^) — dki, 
where probabihty u;^(k, ki) is obtained from u;|(k, ki) with replacement Al"(k, ki) by 

M^(k,kO = ^a^f^T«'*'M[nkjA . (9.12) 

k ki n^ko \ ujuji I i i 

In contrast to preceding case, M.'^ represents the sum of two contributions. The first 
one defined by bare three-gluon vertex is connected with self-interaction of a gauge field 
and has no analogy in Abelian plasma. The second contribution is associated with the 
scattering of wave by the screening cloud of the charges. 

The remaining terms in the r.h.s. of (9.8) and (9.11) represent the interference of 
Thomson scattering by longitudinal and transverse virtual waves, respectively. It is easily 
to see this, having used 

Im(5r^W,ki) = -Snul^ul J ^ ^^^^^ S{vh), 

lm{*r^'k%[nk2f) = -SttuI^u, f ^ M^([nk2]v) <5(t;A;2). 

^ J 47r VK 

The interference of the last two mechanisms of scattering is absent. 
Thus, summing the preceding, instead of (7.1), (7.2) we have 

■wi 



where 



^\k)^?,{uj,i/gf I (4-Og(k,kO(-^)dki, (9.13) 



I —S{u;l - < - v(k - ki)) I A^'=(k, ki) + A^ll(k, ki) + M^{k, k,: 
The interference term in (9.14) between At" and M.-^ vanishes by relation 

ki[nk2y f ^vVS{vk2) = 0. 



It is convenient to interpret the terms intering to = M.'^ + A4" + M. using a 
quantum language. In this case the term Ai'^ connected with the Thomson scattering can 
be represented as the Compton scattering of the quantum of the soft modes (plasmon) 
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by QGP termal particle. defines the scattering of a quantum oscillation through 

a longitudinal virtual wave with the propagator A\k2), where a vertex of a three- wave 
interaction is induced by HTL-correction ST 3. A4-^ defines the quantum oscillation scat- 
tering by transverse virtual wave with propagator A*(A;2). In this case A4 fulfils role of 
the total scattering amplitude. 

At the end of this Section we note the principal distinction of obtained expression of 
the nonlinear Landau damping rate from the damping rate for hard particles [31]. In the 
last case the scattering amplitude involves the resummed gluon propagator in the electric 
and the magnetic channels only, i.e. A\k2) and A*(/c2) respectively. In our case, for small 
particle momenta (|k| < gT), the scattering amplitude takes into account the thermal 
masses of the particles and exchange contributions (vertex corrections). This leads to 
more complicated expression for A4 involving the effective vertecies, in contrast to similar 
expression for the fast particles. 

10. THE ASSOCIATION WITH HTL-APPROXIMATION 

The kernel Q{k, ki) possesses two main properties. The following inequality results 
from definition (9.14) 

g(k,ki)>o. (10.1) 

Next from (9.4) it follows that A1'^(k, ki) = 7Vf*'^(ki, k). The correctness of this equality 
results from the conservation law (6.13). For A^" and TW"*" we have the similar relations: 
Al"'^(k, ki) = (All''-'-(ki, k))*. Their proof trivially follows from the definitions of 
and A''*. The consequence of these equalities is a main property of a symmetry of kernel 
(5(k, ki) with respect to permutation of a wave vectors k and ki 

Q(k,ki) = g(ki,k). (10.2) 

Prom (9.13) follows that in the case of a global equilibrium plasma and by inequality 
(10.1), waves of a high frequencies are damped out, and a smaller ones are increased. In 
particular, in the limit of |k| ^ we obtain from (9.13) 

y(0) = 3{u;,i/grj (cu^i - u;lJQ{0,k,) ^-^j dk^ < 0, 

i.e. k = 0-mode is not damped, as it was calculated in Ref. [17]. 

This clearly shows that nonlinear Landau damping rate (9.13) which is not of fixed 
sign, cannot be identified with gluon damping rate, calculated on the basis of resumming 
Braaten-Pisarski techniques. Actually 7'(k) defines two processes: the effective pumping 
of energy across the spectrum sideways of small wave numbers with conservation of exci- 
tation energy and properly nonlinear dissipation (damping) of longitudinal plasma waves 
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by QGP particles, where the first process is cruciaL To see this, we rewrite the kinetic 
equation (6.14) as follows 

^^-?>{u,,/gf j (4-O4Q(k,k0 (^^J [^jdki. (10.3) 

Hereafter, we have restricted ourseves to spatially-homogeneous case. 

We perform replacement o;^ = {oj\^ — <^k;J/2 + {oj\^ + a;[.J/2, then last equation can be 
presented as 



dWl 



^-J Qs{k,k,)Wiwi^dk, - J QA{kM)Wiwi^dk,, (10.4) 



dt 

where we introduce symmetric and antisymmetric kernel, respectively 



QsikM) ^ 3(^,,/^)2^4-rT^Q(k,ki), 



Qa{K ki) = Hu;,i/gf M^^_K)!Q(k, k^)- (10.5) 

We integrate Eq. (10.4) over dk and introduce the total energy of longitudinal QGP 
excitations: Wq = J W^dk. Then by the properties of kernels (10.5) we obtain 



dW}, 



J Qsik,k^)WiWi^dkdk^<0, 



dt 

i.e. properly nonlinear dissipation of excitations is defined by symmetric part of a kernel 
Qs only. It is necessary to compare the parts of nonlinear Landau damping rate which 
correspond to the nonlinear dissipation of waves only, with damping rate of boson modes 
from HTL-approximation, namely 

7^(k) ^ 3Mgr f ^^TTX^ ^(k, kO<dki. (10.6) 

The function 75(k) is positive for any value of a wave vector k and particular, for |k| = 0. 

Antisymmetric part of a kernel is not associated with dissipative phenomenon and 
defines the spectral pumping from short to long waves. It is easy to see, by considering 
the model problem of interaction of two infinitely narrow packets with typical wave vectors 
ki,k2. 

Let us introduce as follows 

Wiit) = Wi{t)6(k - ki) + W2{t)6{k - ks), |ki| > |k2|. 
Substituting the last expression into (10.4), we obtan the coupled nonlinear equations 

= -{Qs + Qa)WiW2 , W,{to) = W,o, 
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-gf = -{QS - Qa)W,W2 , 1^2(^0) = W2O. 

Here, Qs,a = Qs,a0^1: ^2)- The general solution of this system has the form 

g-C(i-io) 

W,(t) = -W,o C_^^^ ^ ^^^^^^ ^ _ ^^^^^^ , 
W2{t) = W,{t) - ^ , C^{Qa- Qs)W,o + {Qa + g5)W-2o. 

^5 + Qa Qs + Qa 

For |ki| > |k2|, by the definitions (10.5), C > and therefore, in the hmit for i — > 00 we 
have 

V5 + ^A 

Here Wq — Wiq + W2Q is total initial energy of packets. Thus we see that as result of 
nonlinear interaction of two infinitely narrow packets the effective pumping of energy 
across the spectrum sideways of small wave numbers takes place. 

In this result of the pumping, the part of excitations energy (proportionaled 

to Qs) nonlinear absorbs by QGP particles. The absorption value converges to 

AM^=«-0(W^io/0. 

The process of nonlinear scattering of plasmons by QGP particles only, not results in 
their relaxation in homogeneous isotropic plasma. In fact, by the kinetic equation (10.3) 
the general plasmons numbers conserves 



dN^ d 



dt 



d r 
" diJ 



Nldk = 0. 



It follows that if for time t — to, Nq plasmons are excited in a plasma, then excitation 
energy for any t > to does not less than a value oOpiN^ by conservation law of plasmon 
numbers. In homogeneous plasma the total dissipation of longitudinal excitations energy 
is defind by slow processes of four-wave interaction. 



11. THE ESTIMATION OF 7^(0) 



Now we present a complete calculation of 75 (k) at zero momentum of an incident 
field. We start from representation of 75(k) in the form of (10.6) with kernel (9.14). We 
introduce the coordinate system in which axis OZ is aligned with vector ki; then the 
coordinates of vectors k and v equal k = (|k|, a, /?), v = (1, 6, (f), respectively. By $ we 
denote the angle between v and k: (vk) = |k| cos$. The angle $ can be expressed as 

cos$ = sin 6* sin a COS ((^ — /3) + cos 6* cos a. (H-l) 
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In the limit of |k| — > 0, the kernel (9.14) reduces to 

g(0, ki) = T^N^ujpi^-^^ j — 5{cos e - pLJ l-M'(0, ki) + M\\ (0, ki) + A<^(0, ki) 1^, 

(11.2) 

where 

< = </|ki|, kl = «)2 _ pi^^ = (4^ -u;pO/|ki| > 0, dO = sinMM(^. 

By using the definitions (9.4), (9.9) and (9.12), we find expressions ^^'^(O, ki), Al"(0, ki) 
and Ai^{0, kl). The first of them is defined more simply. In the limit of zero momentum 

Calculation of A^ll(0, ki) is more complicated. From (9.9) we obtain 

M^OM) = (—] lini -L(^A\k2)6r'\k,-ku-k2)k'kik'^). (11.3) 

\ujpij ujI^ kipl^ |k|-o|k| 

Using the definition of HTL-correction (7.6) to bare three-gluon vertex, after shghtly 
cumbersome computations we define 

|lm ^5r*^''(^, -'^i, -k2)k'k[k{ = (11.4) 



^, k[ {pj.f ^^^^^ 

'ujpi 1- (pLJ2 



- cos ct— — ^^yif-^ lim A ^'(/c2) + Slkil^cosapk^Vkj. 



In the last equility we use the definition of A~^'(/c) as the function F(a;/|k|) (6.9). In- 
serting (11.4) into (11.3), we reduce the scattering amplitude with longitudinal virtual 
oscillation to 

M\QM)^-{—\ ^i^cosa + 3f^) ujpi{l-{pi^f)cosa\im A\k2). (11.5) 

Now we consider the hmit of the term with transverse virtual oscillation. Prom (9.12) 
it follows 

M\^M) =f^V^ ^"^ hm -^(A*(fc,)T-^;t^F,[nk,]0. (11.6) 

\^pij <^ki kfsmo; |k|-^o|k||n| 

Using the definition of effective three-gluon vertex (7.4) and the relation 

6U}pi |k|^0 6UJpi 
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we obtain 

lim *r^\k, -ki, -k2)k'ki\nk2Y = 

|k;|->o |k||n| 

Q 1,4 1.2 I 

= -- -isina (1 - - ^ ,lrn^""^^2) ^11-^) 

Substitute (11-7) into (11-6) and take into account the relation (11.1), we obtained 
required limit 

M^{Q,ki) = - [—] ^sin^sinacos((^-/3)- (11.8) 



o \ , I |ki|(«) '-pLjsin^sinacos((^-/3) lim A*(/C2). 



A (a. 

2 \ujpi 

The terms in the amplitude A4 not containing A''* combine to give 

( — ] {cos^cos^ — (pkj^coso; — p^^ sin ^ sin a cos ((^ — /?)}. 

By the 5-function in (11.2), the expression in curly braces vanishes. Thus, all terms in 
Ai, not containing the factors A''*, arc in relatively reduced in the limit of k = 0-mode. 

The remaining terms, after substitution into (11.2) and integration over solid angle, 
yield 

q AT 1.2(1 _ (J \2\ 

g(o,kO = W-f )\ ^^1^^ > vl^e{i - pLJ{(i - (pLJ^) cos^^i lim (a'(^,))|V 

+^ii<r'-pU" sm'a\ hrn (A\k2)m. (11.9) 

We note that this expression is not dependent on angle /3. This enables us to represent 
the integration measure in the r.h.s. of (10.6) in the form 



OO —1 

J dki ^"271 J k^d|ki| J d(cosQ). 



Futhcr we suppose that the excitations in QGP become isotropic over the directions of 
vector ki on a time scale which is much less than time scale of the nonlinear interaction. 
This enables us to consider the spectral density M^^i ^ ^ function of |ki| only. Now we 
introduce the spectral function 
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such that the integral /q°° W^|',^^|d|ki| = W'' is total energy of longitudinal oscillations in 
QGP. Substituting (11.9) into (10.6) (for k = 0) and performing the angular integration 
over a, we obtain finally 7^(0) 

\K\ 

7^0)- J g(|ki|)iy|L,|d|ki|, (11.10) 



where a kernel (5(|ki|) has the form 

g(|k,|) = ^Tr/iV, ^)^^ ; ^^^^P IkiKpLj^ ^(1 - pLJ{(1 - (pLJ') I hm {A\h))\'+ 

-pLJ'I ^^rn m^m, (11.11) 

and the upper cutoff \k.l\ distinguishes between soft and hard momenta: gT <^ |k^| <^T. 
For crude estimation of 7^(0) let us expand the kernel (5(|ki|) in the momentum |ki|, 
using the approximations 

I I 3 ki I ujpi 

"-^""^"^-^"lo;:^;'"-^"^' ■■■ ■ 

Keeping the leading in |ki| term in the expansion (11.11), we obtain 



The function VT^i is approximated by its equilibrium value [17]: PF^i ~ ^ttT and therefore 

W(^^\ ^ len^kjT. (11.13) 

Substituting (11.12) and (11.13) into (11.10), and setting |k^| ~ gT, we finally define 

7^(0) ^ +1.047Ve/r. (11.14) 

In our case the coefficient of the g^T has the same sign but is significantly large than 
corresponding one, calculated in [4]. 

12. THE 7'(0) DEPENDENCE ON A GAUGE PARAMETER 

Now we estimate the contribution of a gauge dependence term of a propagator in 
covariant gauge (6.8) to the nonlinear Landau damping rate 

'^Uk2) = ^k2pk2a/[{uJ2 + i^f - k^]^, = k - ki. (12.1) 
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By the Ward identities (7.10), a gauge-dependent part of 7'(k) can be represented as 



'""'g,^;';"'' (A-■'W-A-^'(^.)) ^, , (12.2) 

If we take into account that A~^''(k)\^^^i = A~^'(A;i)L = 0, then formally, a gauge- 

dependent part of nonlinear Landau damping rate vanishes, as it was mentioned in Sec. 8. 
However, we show that the integral in the r.h.s. of (12.2) develops on mass-shell poles. We 
consider, for example, the coefficient of (A~^'(/c))^. In the limit |k| — > this coefficient is 
equal to 

J V du, ^^=,.^''^[(a;^,^-a;,,-ie)2-k?]2 

(numerical factor is omitted). For lower limit, integrand expression (for Wl^^ ~ const) is 
as follows 

d|ki| 

|ki|~0 

i.e. the integral involves a power infrared divergence and is infinite at the pole. Thus, 
a gauge parameter in (12.2) is multiplied by x oo uncertainty. We investigate this 
uncertainty, following by reasoning [27]. 

In order to evalute the double poles in (12.1) we use the prescription 

1 ,.5 1 

lim 



[{uj2 + iey - ki]2 m2^o dm^ [{lu2 + ie^ - ki - m^] 

d 1 / 1 1 \ , , 

= hm , , , . 12.3 

m2^o dm' 2^\4 + w? ^UJ2 - Vki + m2 + ie uj2 + \J\^l + + ie^ 

By computing we first take befor the effective on mass-shell limit uj Upi. 

We focus on the damping rate 7^ of excitation at rest (at vanishing three- momentum). 

Performing the interesting limit |k| — > in (12.2) and taking into account (12.3), we 
define 

V(o) = ^e/iVe- hm ^ J^^(Wi^\(dR^y^ (±y 

^ 3 "upi w?^Q dw? J ^kf+m^ duji J uji=ujl^ \ kl J 

(A-^^(a;, 0) - A-i'«, kO)'{<5(a; - < - ^kf + m?) - 5{uj - < + .Jn + m?)]. (12.4) 
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Without restriction for the general case we choose u > Upi + m. Then we obtain from 
(12.4) 

^« 3 ^ ^Upi m^o dm? \\^klQ + V ^fcio ^ ^ ^^i / '^i=<„ [Kio)^ - ^lo]^ 

(12.5) 



(A-i'(a;, 0) - A-i'(a; - ^kl, + m\ + 
where fcio is a solution of equation 



dkl akfy|ki|=fcio 



< = - V + ^ ■ 

In order to perform the interesting limit w? and uj — > a;^;, we notice that the solution 
kiQ vanishes for w? = as 

fcio c^u -ujpi + Oikl^). 

Working out the derivative with respect to 8/ dm? and taking w? — > 0, we find that the 
most singular term in the limit kiQ comes from the derivative 

a / fcio \ 1 



Using the approximation Wl_^^ a; AnT and A~^'(a;, 0) — A~^'(a; — /cio, /cio) ~ 2a;p;A;io, we 
finally obtain 

iiO)c.-2^^^g^N^T[^). (12.6) 

From the last expression we notice that by going on mass-shell, the gauge- dependent part 
of the nonlinear Landau damping rate vanishes. 

Excess factor (a; — ujpi) of numerator (12.6) is arised from the function 

dkl 9k?y|ki|=fcio 

in the expression (12.5). In paper [27] this factor is compensated by singularity of statistic 
factor 

T 



lim lim fg{J kiQ + m^) 



associated with the spectral representation of bare propagator l/{k^ + m^). In our case 
this factor is absent. 

One can lead to the result (12.6) in a different way. In derivation of the expression for 
7'(k) (6.15) we use representation of the spectral density I^^ in the form (6.10). Gen- 
eral speaking, this representation holds in the linear approximation only, when the time 
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correlation (dependence on uj) is one-to-one correspondence with the spatial correlation 
(dependence on k). To include the effects of wealcly nonstationary ingomogeneous plasma 
motions (when one-to-one correspondence between excitations frequency and its wave 
number fails), we replace the sharp 5-function in the Ij^y^ by Breit-Wigner form [6, 32, 
29], with width V 

ji =—(p 1. \- p 1. ^ (^9 7) 



The parameter 7^ is interpreted as damping rate of boson mode of order g^T. Note that in 
this case we are dealing with correlator {A'^^{Xi)A^j^{X2)) dependence on difference X1—X2. 
Dependence on the midpoint {Xi -\- X2)/2 is accounted in the form of a dependence on a 
slow coordinate and a slow time in the l.h.s. of kinetic equation (6.14) by means of the 
operator d/dt -\- Y{p/d-x.. 

Using representation (12.7) and formula 

1 (91/1 1 \ 



[{uji — ujpi — ieY — kf]2 9kf 2|ki| \uji — tOpi — |ki| — ie uji — uipi + |ki| 
in this case we derive, instead of (12.5) 



It 



= -^;/iv,-3-/ k;d|k.i / d.. _ ^J.^, ^ 

(A-"(a.,0) - A-"(a.i, 0)f{6(u,i - a.,, - |ki|) - 6(u,i - + |k,|)}. 

Perform the integration with respect to over da;i, and differentiation with respect to |kip, 
we obtain expression 



,.(0).-4^ifL,,.,.T(^)/k?d|k. 
It follows the result (12.6). 



13. CONCLUSION 



Let us consider in more detail approximations scheme, which we use in this paper. 
In fact, here two types of the approximations are used. The first of them is connected 
with the employment of usual approach, developed in Abelian plasma, to QGP, i.e. the 
standard expansion of the current in powers of the oscillations amplitude and computation 
of interacting field in the form of a series of a perturbation theory in powers of a free field 
(more precisely, in gA^^^). However, in contrast to Abelian plasma, in our case even if 
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the first two nonlinear orders of the color current are taken into account, much more terms, 
defining the nonlinear scattering of waves are derived. Here, we use second approximation, 
connected with the notions going from the papers devoted to high-temperature QCD or 
more precisely, the set of rules for power counting these terms, developed by Blaizot and 
lancu [12]. These rules enable us to singled out the leading terms in the coupling constant. 
These terms are purely non-Abelian in a full accordance with the conclusions of Ref. [12]. 

However we note that in [12] behaviour of a mean field (A^(X)) is investigated, which 
in our case vanishes. In the covariant derivative "D^, containing only the random part 
of a gauge field in our approach, we suppose that d/dX^ and gA^ are of the different 
orders. In field strength tensor we distinguish the linear and nonlinear parts. Therefore 
the aproximatons scheme which have been carried out in this paper suffers from disad- 
vantage of breaking non-Abelian gauge symmetry of a theory at each step of approximate 
calculation, as it was discussed in Sec. 4. 

This is really true if we suppose that the magnitude of a soft random oscillations 
|74^(X)| is of order T or |A^(A;)| ~ l/g{gTY in the Fourier representation. In this case, 
by using obtained expressions for terms in the expansion of a colour current (3.2), (5.4), 
(5.9) and estimations (5.7), (5.12), we have 

i.e. all terms in the expansion (2.14) are of the same order in magnitude and the problem 
of resummation of all the relevant contributions appears. Thus a gauge symmetry is 
recovered. 

In this paper we have restricted ourselves to just a finite number of terms in the 
expansion (2.14). This impose more rigorous restriction on the magnitude of random 
oscillations: \Al{X)\ ~ gT {\Al{k)\ ~ l/{gTf). In this case we have 

each following term in the random current expansion is suppressed by more powers of g 
and use of the perturbation theory is justified. A gauge symmetry is restored when we 
take into account consistently all contributions at the leading order in g (Sec. 6) to the 
probability of the soft excitations scattering by the hard thermal particles. In this case we 
derive the expression for the nonlinear Landau damping rate (7.1), (7.2) which is closely 
allied in the form to corresponding gluon damping rate in HTL-approximation [4]. 

In Sec. 10 it was shown that the nonlinear interaction of longitudinal eigenwaves 

leads to effective pumping of energy across the spectrum sideways of small wave numbers. 
Consequence of this fact is the inequality 7^(0) < 0, i.e. k = - mode is increased. 
The dissipation of the energy of plasma waves by QGP particles here, is not lead to 
total relaxation of a plasma excitations. In the scale of a small |k|, effects described by 
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nonlinear terms in the expansion of the colour current of higher-order in the field, come 
into play. Consideration these effects involves suppression of increase of k = 0-mode. 

The r.h.s. of obtained kinetic equation (6.16) in the regime |k| ^ gT contains the 
part of possible processes of a plasmon scattering in QGP only, namely, the processes of 
a type 

9* + 9^ 9* + 9 , 
9* + ^9* + , 

where g* is plasmon collective excitations and g,q,q are excitations with characteristic 
momenta of order T. Diagrammatically this corresponds to graph with four external lines, 
where one of the incoming (outcoming) lines is a soft and the other is a hard. However 
it is clear that there arc further contributions to the damping rate of a soft gluon, going 
without exchange of energy between a hard particles and waves. They are associated with 
the processes of nonhnear plasmon scattering by a soft excitations of QGP, i.e. with the 
processes of a type 

9* + 9* ^9* + 9*, (13.1) 

9* + q*{q*)^9* + q*{r), (13.2) 

where q*, q* are plasmino collective excitations. The kinetic equation describing the pro- 
cess (13.1) and (13.2) is the equation of a purely Boltzmann type, i.e. collision term in the 
r.h.s. of this equation has standard Boltzmann structure, with on gain term and a loss 
term. The probability of these processes is defined by preceding methods from nonlinear 
current of a fourth order jj^^'', if the process of interaction iteration of higher-order in the 
field is taken into account. 

As it was mentioned in Introduction, the Boltzmann equation was already used for 
definition of damping rate of the fast particles. In paper [11] (see also [14]) on the basis 
of Boltzmann equation, the damping rate for hard gluons in the leading logarithmic or- 
der has been computed. The value of obtained the damping rate is fully coincident with 
corresponding damping rate derived in quantum theory [6, 7]. The scattering matrix 
element appearing in collision term corresponds to elastic scattering of hard gluons in the 
resummed Born approximation. By Heiselberg and Pethick [11] was noted that the par- 
ticle damping rate is considerably more difficult to calculate for small particle momenta, 
|k| < gT. Particles momentum becomes of the same order with momentum transfers. As 
in the case of scatterng of plasmons by hard particles discussed above, vertex corrections 
should be taken into account also. It leads to more complicated expressions for probabil- 
ities of both plasmon-plasmon (13.1) and plasmon-plasmino (13.2) scattering in contrast 
to (9.14). 

A detailed research of the process of (13.1) and its infiuence to relaxation of soft Bose 
excitations in QGP will be presented somewhere [26]. Here, we note only, that it is the 
process of higher-order in the field, which probability is defined with the help of three- 
gluon, four-gluon effective vertecies and effective propagator, as the probability of the 
process of nonlinear scattering of a plasmon by QGP particles, derived in this paper. 
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APPENDIX A 

We rewrite the system of equations (9.6), (9.7) in the coordinate represetation. In the 
temporal gauge we have 

d^x g 



^VF^Q"(E" - v(vE")), 



dt^ m 

yr^^(vA'')Q^ (^1) 



dt 

where i is a coordinate time; v = dx/dt, E" = —dA°-/dt, and 

A"(x,i) = j A^e^^'-^^^k^dk , A^ = i^A^. (A2) 

The Eq. (9.6) corresponding to the component A* = 0, becomes identity. 

If we neglect by the wave action on typical particle, then the colour charge motion will 
be constant and its colored vector will be fixed with initial condition: Xo(t) = Vot, Q" = 
Qq. To derive the oscillations of a color particle, excited by fields and being linear order 
in amplitude of field, it is necessary to neglect by variations of v, Q" in the r.h.s. of (Al) 
and set 

A"(x,i) ~ A"(voi,i) = j A^e-^^'^k-vok)*^!^, 

instead of {A2). 

In this approximation the solution of system {Al) has the form 

x(t) = vot + {-^) QS0^ / 7 "'^"t? e-^(-l^--^)* dk ^ xo(t) + Ax(t), 

{t) = + igr'^Ctl j e-^('^^-"°'^)* dk = + AQ" {t) . (^3) 

J uj]^ — Vok 

Using these expressions we derive the radiation intensity of oscillating colour charge. 
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It is equal to work of radiation field with charge in unit time 

= / (E^jQ)dx, {A A) 

where Eq is a field induced by a colour current jg of a charge Q". The sign in the r.h.s. of 
(A4) corresponds to choose of a sign in front of current in the Yang- Mills equation (2.2). 

We introduce Eq and Jq in the following form 

E^(x,i) = j E^^^e^-^^-^^kdo;, j^(x,i) = j j^^^e^"^^--* dkdo;. {Ah) 

The Fourier-component of a field E^^ = k£^k,a;/|k| is associated with j^^ by Yang- Mills 
equation 

Substituting (^45) and (^46) into (A4) and taking into account reahty of W', we obtain 

= ¥ / ' (zdm - ro) ^«si^;.>e-<"'-""d.d.'dk. (^7) 

Here, by random character of a current phase, we replace the product Jk"c!;ik^a;' by its 
averaged value. 

The colour current formed by the constrained motion charge {A2)) is presented as 

dk 



j^(x,t) = 5v(t)Q'^(t)<5(x-x(t))=5v(t)Q'^(t) / e'i'(---W)^ = / j^(t)e'^Mk, 
where we have in the linear approximation 



m = -i-Mm\t) e-'^^w ~ (^8) 



(27r) 



and Ax(i), Av(i), AQ" are defined by (A3). 

The first term in the r.h.s. of (^48) is connected with hnear Landau damping which is 
absent in QGP. Because of this fact this term is omitted. The term in braces, proportional 
to Qq yields the Abelian contribution to radiation and is connected with usual spatial 
oscillation of a colour particle. The term with AQ" gives the non- Abelian part of radiation 
and is induced by precession of a colour vector of a particle in the field of incident wave. 
As it is easily to see the interference of these two contribution vanishes. Further we shall 
restrict our consideration to the second part of radiation which induced by non- Abelian 
part of a colour current 



i,^'k.,uj) non— Abelian j (jk(^))non— >lfeeZian G 
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is' /-QS / 5{u - ul + Vo(k' - k)) dk'. 



(2Tif ^"i 4, -Vok' 
Substituting the last expression into (A7), taking into account 

and integrating over dcj, dc^;' and dk", we obtain with replacement of argument k' ki 

Here, gr^ = Q^Qo; '^ki = <^ki ~ (^i — k)v and the suffix "0" of velocity is omitted. 

For imaginary part of longitudinal permittivity we use approximation similar to [21] 

Going from the spectral density /^i to density of longitudinal oscillations number A^^i 
setting the constant (f equal to 2/(27r)^, we can cast radiaton intensity expression (A9) 
into the final form 

r dkdki 

{y\^%on-AbeUan = j <(k, kO^^"^, (^10) 

where w;^(k, ki) is defined by (9.3). Expression (AlO) represents radiation intensity of 
isolated colour charge moving in a quark-gluon plasma in the direction to v of a field of 
a longitudinal wave with frequency uj\^ and wave vector k. 
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